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Abstract: We study the effects of strong magnetic fields (B > fO 13 G) in the cross 
section for v e n — > pe scattering in the presence of a degenerate electron background. This 
can be relevant for the v e propagation in the proto-neutron star stage after supernovae 
collapse. We find that for field strengths B > 10 16 G(E U /10 MeV) 2 the v e opacity is 
sizeably affected by the magnetic field and can lead to a shift in the location of the 
electron neutrino sphere towards lower densities. We discuss the implications that this 
may have for scenarios proposed to explain the observed pulsar velocities. 
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There has recently been a renewed interest in the study of the neutrino propagation 
at very high densities (p > 10 10 g/cm 3 ) and in the presence of strong magnetic fields 
(B > 10 13 G), in connection with the proposed explanations for the observed pulsar 
recoil velocities, which could arise from an asymmetry in the neutrino emission during 
the first seconds after the collapse of massive stars [(j], [| [|, |], H 0> |U • The high densities 
reached in the proto-neutron star stage result in the trapping of the neutrinos within the 
so-called neutrino spheres, i.e. the surfaces from which the optical depth becomes of 
order unity. The location of these surfaces depends on the neutrino flavour because of 
the effects of the charged current interactions, which allow the v e to escape only from 
larger radii (mainly because of the large v e n — > pe cross section) than the v^ T (and their 
antineutrinos), while the u e -sphere radius is in between the previous ones (due mainly to 
the v e p — > ne + interaction). 

The fact that the radii of these neutrino spheres is flavour dependent is at the basis 
of the mechanism suggested by Kusenko and Segre || to account for the pulsar kicks, 
which is based on the resonant conversion v e — > v T just in the region between those two 
spheres, so that trapped v e s oscillating to z/ r s become free to escape. As was noticed by 
those authors, the large magnetic fields which may be present in the proto-neutron star 
(as evidenced by the large B fields, B ~ 10 12 -2 x 10 13 G, observed in young pulsars), may 
be responsible for a distortion of the resonance density for the MSW conversion and can 
lead then to an asymmetry in the average energy of the emitted v T s. For this mechanism 
to be successful, field strengths B > 10 14 G would be required, and some scenarios have 
been proposed in which such large fields may be achievable at the epoch of the supernova 
collapse [y, [10], [11], [12]], with the maximum conceivable strength being B ~ 10 18 G, a value 
beyond which the magnetic energy becomes larger than the gravitational one. 

We want here to study whether these large magnetic fields may actually affect the 
v e n — > pe cross section which determines the v e opacity itself and hence shift the location 
of the z/ e -sphere. This can happen due to the modification of the phase space distribution 
of the final state electrons in the presence of the magnetic field. It has indeed been shown 
that the neutron decay and electron capture processes |TJ, [14|, |l5j are affected significantly 
for B > B c = m 2 e /e = 4.141 x 10 13 G, in which case the spacing between the first Landau 
levels becomes larger than the electron mass m e . For larger B fields, the separation 
between Landau levels may become comparable to the typical energies of the neutrinos 
emitted in supernovae, E Ue ~ 1-30 MeV, and hence the neutrino opacity may also be 
changed. 

The reaction rates for the v e n interaction in the presence of magnetic fields were 
considered in [plj, focussing in the conditions for (3 equilibrium. Furthermore, the possible 
B dependence of the v e n — *> pe cross section has been recently invoked as an explanation 
of the pulsar kicks if an asymmetry in the magnetic field distribution develops during the 
first second after the collapse [I7|. A detailed evaluation of the magnetic field effects on 
the v e opacity is however lacking and will be the subject of the present work. 

We will then concentrate in the study of the u e produced during the Kelvin Helmholtz 
cooling phase |TS|, [T9| , just after the time when the shock produced by the core bounce 
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dissociates the heavy elements initially present and make the v e opacity in the shock 
heated region to be dominated by the reaction v e n — > pe (before that, scattering off 
nuclei dominates). 

At zero density and in the absence of magnetic fields, the total cross section for this 
reaction is (see e.g. [|(J ) 

G 



°o = —{3v + Wa){E v + QYJl - m\j{E v + Qf 



7T 



where G = Gfcos6 c stands for the product of the Fermi coupling and the Cabibbo 
mixing while gy — 1 and qa = 1-26 are the vector and axial nucleon couplings. Since the 
energy transfered to the recoiling proton is 0(E^/m p ), which is negligible for the neutrino 
energies we are interested in, one has E e ~ E v + Q, with Q = m p — m n = 1.293 MeV. 
This 'elastic' approximation is indeed quite reliable pl| for densities below the nuclear 
saturation one and for non-degenerate neutrinos (fi u — 0), as is the case of interest here 
for the study of the location of the neutrino-sphere. 

At finite densities (but still with B — 0), the main modification will come from an 
overall electron blocking factor 1 — f e {E e ) multiplying the cross section, where 

fe{E) = i + expp;- M )/r]' (2) 

is the Fermi Dirac distribution at temperature T, with /i being the electron's chemical 
potential. 

Considering now a non-vanishing magnetic field, the matrix element for the process 
will remain essentially unaffected (we average over initial spins neglecting the neutron 
polarisation, which is small for B < 10 17 G) [0, The main modification will then 



come from the available phase space for the electrons, since the phase space factor for 
B = 

d 3 p 



has now to be replaced with 

? = (£F (4) 

where go — 1 and g n = 2 for n > 1 are the degeneracies of the Landau levels of energy 



E n = ^/pl + ml(l + 2nB*), (5) 

and we have introduced £>* = B/B c . 

Computing the v e n pe cross section in the presence of the magnetic field we then 
obtain 



* = ib(9v + ZglXB* E n 9nJdp z 6(E u + Q - E e ) [1 - f e {E e )) 

= i(9v + Z9\WeB* [1 - fe(E„ + Q)\ (E„ + Q) En9n/P^, (6) 
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where pf = y \E V + Q) 2 - m\{\ + 2nB Sf ). The maximum Landau level accessible to the 
final state electron, for a given energy of the initial neutrino, is 



n, 



int 



25* 



E u + Q 



(7) 



In figure 1 we plot the cross section in vacuum, normalised to the B = one (E 
c/cq), as a function of the magnetic field, fixing E u = 10 MeV for definiteness. 




Figure 1: Normalised cross section, £ = ct/ctq, for i£„ = 10 MeV, as a function of the magnetic 
field = B/B c , in the absence of background matter and for T = 0. 



For field values 



B > 



E v + Q 



B r ~ 10 16 G 



E v 



10 MeV 
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only the lowest Landau level (n = 0) contributes to the phase space, and in this case one 
has 2 

+ 1g^m\B^ (9) 

which grows linearly with B and is independent of the neutrino energy. For smaller 
magnetic field values, n max > 1 and hence more Landau levels contribute to the sum 
in eq. (||). The singular behaviour present each time that a new Landau level opens 
up, which is similar to the one found in the /3-spectrum in n-decay []14|| , arises from 
the p^ 1 factor from dp z = E/p z dE, and is expected to be somewhat smeared once the 
proton recoil momentum is included and its effects also averaged out once a distribution 
of neutrino energies is considered. In the limit of small magnetic fields, n max will be large 
and we may approximate the sum over Landau levels as 



"•max 

n=0 Pz 



{E 2 e - ml)- 1 / 2 + 2YZ=T( E e ~ + 2nB *)) 



-1/2 
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i + Kfk /x dy(l-y)- 1/2 



~ {E 2 e -ml)- 1 ' 2 

where we defined x = (E e /m e ) 2 — 1. From this we get, for Z?* <C x, 

„ -B* 

E ~ 1 -, 

2x 

and the 5 = result in eq. (Jj) is then asymptotically recovered. 

In the presence of background matter, the electron density is given by 



(10) 



T,fe(E) = ^Y,9nJd Pz f e (E n 



(2vr 
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and a similar expression holds for the positron density, with the replacement /i — > — /i. 

Taking into account that for typical proto-neutron star temperatures a positron back- 
ground is also generally present, we can introduce the electron fraction as 



Y, 



The star density may then be written as 



n e - n e + 



Tip ~\~ Tl n 



VTir, 



p ~ m p (n p + n n ) = —(n e - n e +). 



(13) 



(14) 
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Figure 2: Normalised cross section vs. background matter density, fixing E u = 10 MeV and 
for mag netic fields = 10 2 , 10 3 and 10 4 . For < 10 2 the results are almost insensitive to B. 

Considering now the effects of electron degeneracy in the v e n scattering, we plot in 
figure 2 the normalised cross section as a function of the density, neglecting for the moment 
temperature effects (T <^ E v ,fi), for different values of B. Under this approximation 
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f e (E) = 9(p — E v — Q), and hence the effect of the background is just to block the final 
state electrons, resulting in a maximum density p max beyond which no scattering can take 
place and a vanishes. For densities larger than p max only neutral currents will contribute 
to the neutrino opacity 1 . Increasing the size of B will modify the size of a as discussed in 
relation with figure 1, and, due to the B dependence of the electron density, also change 
the maximum density beyond which a vanishes, i.e. p max = p max (B, E v ). One has to 
keep in mind that large magnetic fields will also affect the Y e values corresponding to (3 
equilibrium and hence the details of the proto-neutron star evolution. 



E u = 10 MeV, B* = 1 
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log[p(Y e /0.1)] (g/cm 3 ) 

Figure 3: Temperature dependence of the £ vs. p relation, fixing E u = 10 MeV and £>* = 1. 

Turning now to the effects of finite temperatures, we show in figure 3 the modification 
of the previous picture with the inclusion of thermal distributions, plotting E vs. p for 
different temperatures. The main result is that some of the final states with E e > p 
become occupied at finite temperatures, leading to a partial Pauli blocking of the final 
state electrons, while for E e < p, some electron states below the Fermi energy become 
free at non-zero temperature and the available phase space does not vanish abruptly 
but instead diminishes smoothly for increasing densities. It is clear from figure 3 that 
only the low energy part of the neutrino spectrum is affected by thermal effects, which 
just open up some phase space for the final state electrons and result in a non- vanishing 
contribution to the opacity up to much larger densities. 

Let us now discuss the implications of the previous results for the emission of neutrinos. 
During the cooling phase the neutrinos will diffuse out from the inner regions of the proto- 
neutron star 2 and will be emitted from the surface of the neutrino sphere. The presence 

1 Here it is useful to recall that for B = 0, T = 0, one has fi ~ (37r 2 n e ) 1 / 3 ~ 
11 MeV( / 9r e /10 10 g/cm 3 ) 1 / 3 , and hence in this case p max Y e ~ 10 10 g/cm 3 ((^ + Q)/ll MeV) 3 . 

2 Actually, it is the lcpton number which diffuses, since neutrinos can be captured by neutrons and 
then rccmited by electron captures, as well as being thermally pair produced. 
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of large magnetic fields can affect the electron neutrino opacity and hence shift the radius 
of the neutrino sphere. A change in this radius affects the neutrino luminosity in two 
ways: i) it modifies the area of the emitting surface {L v oc R 2 ) and ii) the temperature 
of the matter in the emission region can be different (the temperature associated to the 
neutrinos is usually obtained from a fit to the resulting neutrino spectrum, but the low 
energy neutrinos come from deeper layers of the star, while those of higher energies are 
trapped up to larger radii, making the neutrino spectrum not exactly a Fermi Dirac one). 

Usually a mean neutrino sphere radius is obtained by averaging over the neutrino 
energy distribution, but we prefer instead to work here with an energy dependent neutrino 
sphere radius, since for instance when considering the Kusenko and Segre mechanism one 
needs to actually compare this radius with the location of the resonance for neutrino 
conversion, which is energy dependent. 

An estimate of the matter density at the z/ e -sphere radius R e can be obtained from 
the condition of having unit optical depth, r(i? e ) = 1, with 

/*oo /*oo 

r(r) = / dr ^2niCr(is e i — ► X) ~ a / dr n n (r), (15) 

Jr • Jr 

and where we used that the sum over possible scatterers % is dominated by the scattering 
v e n —> pe. 

Hence, we can write the approximate relation 

l = cr / drn n (r) ~ an n (R e )h n , (16) 

J Re 

where h n = \dhm n /dr\^ is the scale height of the neutron distribution at the z/ e -sphere 
radius (h n is actually not really a constant, since the density profile is not exponential). 
From this one can estimate the density at the neutrino-sphere as 

p(B „)^^ xl0 n^_(H|!X)Vi^)^ (17) 
Y n ah n cm J V E v J \ h n j a 

where Y n = 1 — Y e is the neutron fraction (Y e ~ 10 _1 and h n ~ 10 km typically). 

Hence, a change in a caused by a large magnetic field will shift the density of the 
z/ e -sphere. In particular, for very large fields > 240(^/10 MeV) 2 , for which only the 
lowest Landau level contributes to cr, one has E = (B*/2)(m e /E) 2 , and hence the density 
of the neutrino sphere will behave as B' 1 for large fields. 

Our results differ significantly from those used in ref. [I7|, where the estimate S ~ 



0.775* (for B* ^> 1) was adopted from a simple analogy with the n-decay results |y| . 
As a consequence, to induce a sizeable recoil velocity from the effects on the neutrino 
emission resulting from an asymmetric magnetic field distribution would require much 
larger fields. 

Let us also notice that if p(R e ) is bigger than the p max discussed in connection with 
figure 2 (this can happen for small E u ), the v e n — > pe cross section will vanish, neglecting 
temperature effects, inside the neutrino sphere. In this case, only magnetic fields large 
enough to increase p max beyond p(R e ) may be able to affect sizeably the v e opacity. 
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Considering now the scenario proposed by Kusenko and Segre [EL we note that a 
value of E > 1 would shift the i^-sphere to larger radii, and hence lower densities. This 
may be interesting in order to allow the resonant conversion v e — > v T to happen inside 
the z/ e -sphere but with smaller values of Am 2 (since Am 2 oc N e ), possibly within the 
cosmologically acceptable values corresponding to J2 m vi < 92(Q u h 2 ) eV. This could cure 
one of the main drawbacks of the model, which is the need for large neutrino masses, 
with Am 2 > (100 eV) 2 . However, for typical neutrino energies, E v > 5 MeV, to increase 
S significantly would require quite large fields, B > 10 16 G, and hence this possibility 
seems also difficult to implement. 

Other mechanisms based on large magnetic fields to produce the pulsar velocities, 
such as by means of the B dependence of the differential cross section of URCA processes 
[J], ^ or of ve scattering ||, do not depend directly on the precise location of the neutrino 
spheres, so that the process discussed here will not interfere with those scenarios. 

In conclusion, we have studied the behaviour of the v e n —>■ pe cross section in very 
strong magnetic fields, finding that it can lead to sizeable modifications of the neutrino 
opacities in proto-neutron stars for B > 10 16 G(E U /10 MeV) 2 . The impact of this for 
scenarios proposed to explain the observed pulsar velocities seems then marginal in view 
of the extremely high fields required. 
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